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Enhanced photogalvanic effect in graphene due to Rashba spin-orbit coupling 
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We analyze theoretically optical generation of a spin-polarized charge current (photogalvanic 
effect) and spin polarization in graphene with Rashba spin-orbit coupling. An external magnetic field 
is applied in the graphene plane, which plays a crucial role in the mechanism of current generation. 
We predict a highly efficient resonant-like photogalvanic effect in a narrow frequency range which is 
determined by the magnetic field. A relatively less efficient photogalvanic effect appears in a broader 
frequency range, determined by the electron concentration and spin-orbit coupling strength. 

PACS numbers: 72.25.Fe, 78.67.Wj, 81.05.ue 


Introduction: Two-dimensional electron systems with 
spin-orbit (SO) coupling are currently of a broad interest 
due to the coupled charge and spin dynamics, as revealed 
in a variety of spin related transport phenomena M. 
Owing to the SO coupling, the spin dynamics can be gen¬ 
erated, among others, by a low-frequency electric field (5| 
as well as optically by interband electronic transitions [6|- 
@j. Moreover, an external static magnetic field can enable 
the current generation by light absorption, leading to a 
photogalvanic effect [§] . 

Many of the spin related phenomena, including also 
the ones mentioned above, can be observed in two di¬ 
mensional graphene monolayers and other graphene-like 
materials like silicene for instance. The huge inter¬ 
est in graphene is related mainly to its natural two- 
dimensionality, very unusual electronic structure, and 
high electron mobility which ensures its excellent trans¬ 
port properties [loj, [u|. Even though the intrinsic SO 
interaction in free-standing graphene is negligibly small, 
the Rashba spin-orbit coupling can be rather strong 
for graphene deposited on certain heavy-element sub¬ 
strates EMI. Since the electronic band structure of 
graphene is significantly different from that of a sim¬ 
ple two-dimensional electron gas (2DEG), and the spin- 
orbit coupling creates a gap in the electronic spectrum, 
graphene can reveal qualitatively new effects which can 
not be observed in 2DEG in conventional semiconductor 
heterostructures. Additionally, the SO-related phenom¬ 
ena in graphene are also important from the point of view 
of potential applications in all-graphene based spintron- 
ics devices [15| 17]. 

In this letter we predict an enhanced photogalvanic ef¬ 
fect in graphene. To do this we consider the charge and 
spin currents generated in graphene by optical pump¬ 
ing in the infrared photon energy region, and show that 


the optical pumping can be used to generate in graphene 
not only the spin density [18j, but also a spin-polarized 
net current. An external magnetic field applied in the 
graphene plane plays an important role in the mecha¬ 
nism of current generation. We show that the efficiency 
of current generation per absorbed photon can be very 
high at certain conditions. Apart from this, we also show 
that one can create spin density without creating electric 
current, but not vice versa. 

Model: We consider low-energy electronic spectrum of 
graphene in the vicinity of the Dirac points 111]. Addi¬ 
tionally, we include the Rashba SO coupling [19( and the 
Zeeman energy in a weak external in-plane magnetic field 
B. The corresponding Hamiltonian can be then written 
in the form 


H = hv 0 (±T x kx + Tyky ) + \(±T x a y - T y o x ) + 2(b • er),(l) 


where vq ~ 10 8 cm/s is the electron velocity in graphene, 
t x and T y are the Pauli matrices defined in the sublattice 
space, A = gB is the maximum Zeeman splitting, b = 
B /B, while the + and — signs refer to the K and K' 
Dirac points, respectively. Furthermore, g = gLRB with 
9l = 2 being the Lande factor, and A — oj 2 with a 
standing for the Rashba coupling constant [201 ]. 

The electronic spectrum corresponding to Hamiltonian 
0 consists of four energy bands in each valley, E n k, 
where n = 1 to n = 4 is the band index. The correspond¬ 
ing spectrum for B = 0 is shown in Fig0i, and is given by 
the formula E n ,k(B = 0) = =fA± \/A 2 + fi 2 Wg/c 2 . In turn, 
the expectation value of the spin in the absence of mag¬ 
netic field is oriented perpendicularly to the wavevector 
k [ 20 } , similarly as in a semiconductor-based 2DEG with 
Rashba SO coupling. Contrary to the 2DEG, there are, 
however, no eigenstates of Hamiltonian 0 with a defi¬ 
nite eigenvalue of any spin component. This appears due 
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diagonalization of the Hamiltonian m, and is shown in 
FigU> for wavevectors along the axis x and y. Only the 
states corresponding to the bands labeled in FigJTJr with 
the index 2 and 3 are shown there. As one can note, 
for one propagation orientation the electron bands are 
shifted vertically, i.e. to higher (lower) energy, while for 
the second propagation orientation the bands are shifted 
horizontally, i.e. to left (right) from the point k = 0. 
The latter separation in the k-space of the bands 2 and 
3 is crucial for the enhanced photogalvanic effect. For a 
weak Zeeman energy, A <C a, and for electron momenta 
of our interest, k y/aA/hvo, the field-dependent cor¬ 
rection to the electron energy, calculated by perturbation 
theory, has the form 

E n , k (B) - E n , k (B = 0) = sinfl, (3) 

2 v 0 

where again the upper (lower) sign corresponds to the 
bands (1,4) and (2,3), respectively. 

Assume now that the system is subject to electromag¬ 
netic irradiation. Hamiltonian describing interaction of 
electrons in graphene with the external periodic electro¬ 
magnetic field, A (t) = Aoe _lwt , takes the form 

Ha = ^f'-yo (t • A). (4) 


Figure 1: (Color online), (a) Schematic picture of the band 
structure of graphene with Rashba spin-orbit interaction with 
all possible band transitions for a chosen chemical potential, 
indicated by the vertical arrows. Circles with crosses and 
dots inside correspond to the opposite spin orientations in 
the subbands, (b) Dispersion of the low energy states for 
indicated wavevector orientations. The band index is also 
marked at the plots. The assumed magnetic field B || a; is 
equal to 5 T. 


to the specific form of the SO coupling in graphene. The 
corresponding spin components in the absence of mag¬ 
netic field and for the wavevector k = k (cos 9, sin 9) are 


{&x (k))s=o = T— sin#, (2a) 

Vo 
Vb 

(a y (k)) B =o = ±— cos 9, (2b) 

vo 

where Vk = vo x h vok/ -y/A 2 + h 2 VQk 2 is the absolute value 
of the electron velocity in the absence of magnetic field. 
Here the upper and lower signs correspond to the bands 
(1,4) and (2,3), respectively. Note, both spin compo¬ 
nents give n by Eqs. (l2al) and (12bl) vanish in the limit of 
k = 0 [20| corresponding to the mixed rather than pure 
character of the band states in the spin subspace. 

The electronic spectrum presented in FigQJr is signif¬ 
icantly modified by an external in-plane magnetic field. 
Assume this field is oriented along the axis x. The ex¬ 
act electronic spectrum can be then obtained by direct 


As in Eq. m, different signs correspond here to electrons 
within the K and K' valleys. 

The injection rate of a quantity 0, related to the inter¬ 
subband optical transitions, can be calculated by using 
the Fermi’s golden rule, 


0{u) = Y,° r ' 

n,n' 


o r 


2n 


(ui = 


M, 


d 2 k 


(4/nk | Ha | Hfn'k) 


O" 


xS(E nk + hx- E n ' k )f(E nk ) [1 - f(E n , k )\, 


(5) 


where f(E nk ) is the Fermi-Dirac distribution function. 
Since there are two valleys, K and K', one needs to cal¬ 
culate contributions to 0(w) from both of them. The 
quantities of our interest here are: 

Qn^n’ = (g) 

for the light absorption (where 1 is the identity operator), 


Qn^n’ = - (^kl^l^k), (7) 


for the corresponding spin component injection, and 

O n ^ n ' = <tf„/k|Ji|*„'k> - <*nk|/i|*nk) (8) 

for the current injection. Here, /, = eii, where e is the 
electron charge, while v x = ±vqt x and v v = vqt v . The 
injected spin current, in turn, can be calculated as [2lj 

Qn^n' = _ 0F nk | l^k), (9) 
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Figure 3: (Color online) Normalized charge current I y , for a 
low temperature, T = 1 K (T — 10 K in the inset). Solid 
red line is for a = 4 meV and dashed blue line is for a = 13 
meV. Chemical potential fj, = 5 meV, B is along the i-axis 
and A _L B. 


Fermi wavevector, 


yV 2 + 2A/r A . 

hkp =-=F x— sm 0. 

v 0 2v 0 


( 10 ) 


huj (meV) 


which sets the following boundaries for the Fermi surface: 


Figure 2: (Color online) (a) Normalized charge current I y , in 
the case of low temperature, T = 1 K (T = 10 K in the in¬ 
set). Rashba spin-orbit coupling strength is a = 4 meV (solid 
red line) and a = 13 meV (dashed blue line). The chemical 
potential is fi = 5 meV, B is in the plane of graphene and 
along the *-axis, while A || B. (b) Contributions of indicated 
intersubband transitions to the total current presented in (a) 
for a = 4 meV. 


where j" = \<j v ,i}j\ + / 2. Below we concentrate on the 
results for injection of charge current and spin density. 

Results: Using equations for the injection rate one 
can calculate the charge current and spin polarization 
induced by the optical pumping. Let us begin with the 
photogalvanic effect, i.e. charge current generation. Re¬ 
sults for two different polarizations of the electromagnetic 
held A (t) are presented in Figs. [2] and [3] Here the in¬ 
jection efficiency R is defined as R = RJevpIp, where 
Zb = 'Ke 2 Q/hc and Q is the incident photon flux |22l-l24[. 

The transitions start at Huj ss 2/i — a if /i > a and at 
hfjj « 2/r otherwise. In both cases a strong narrow peak 
in the injection efficiency appears at a resonant energy 
hui « 2/i. This peak is remarkably higher for the electro¬ 
magnetic held polarized along the static magnetic held B, 
compare Figs. [2]and[3j To understand origin of this peak 
lat us consider transitions between the subbands marked 
with n = 2 and n = 3. First, we determine the shape 
of the isoenergetical line corresponding to a given Fermi 
energy, fj, A. For the band corresponding to n = 3, 
one obtains from Eq. ([3]) the hrst-order correction to the 


\! /i 2 T 2A/i A R 2 T 2A/i 

-< n,k F y < - 

v 0 2v 0 v 0 


i//i 2 + 2A/i ^ \/V 2 + 2A/i 

vo v 0 


A 

2v 0 

(11a) 

(lib) 


Due to the k-dependent Zeeman term, the Fermi sur¬ 
face becomes considerably deformed and anisotropic, as 
shown in the left panel of Fig|4] The maximum deforma¬ 
tion is independent of the chemical potential and spin- 
orbit coupling. In turn, the resonance line determined by 
^ 3 ,k — E 2 ,k = RjJ is still a circle given by the condition 


hk u 


y/h 2 io 2 /4 + A hco 
vo 


( 12 ) 


A part of the resonance line is inside the occupied region. 
Therefore, we have an interval of the photon energies, 
(huii, friV 2 ), as shown in the right panel of FigU where 
the transitions occur for positive values of k y , while the 
transitions with negative k y (which would compensate 
partly current) are forbidden. As a result, a very efficient 
current injection occurs in this photon energy window, as 
visible in Fig. [2] In the considered regime of /r A, this 
photon energy interval is determined by the conditions 


huii 


tlUJ2 


2/x — A 


Vk F 

vo 


2/z + A^, 

Vo 


(13a) 

(13b) 
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Figure 4: (Color online). Left: Schematic picture of the Fermi 
line (solid) and resonance line (dashed) for the chosen sub¬ 
bands. Optical transitions are possible only at the part of the 
dashed line outside the filled area. Each transition generates 
a current of the order of 2 eVk F , making the generation highly 
efficient. Right: Side view on the intersubband transitions, 
which can occur in the frequency interval oji < uj < u> 2 - 

which results in the peak width given by the formula, 

h(u 2 -u 1 )= 2^Z±|^A. (14) 

/i + A 

With the increase in temperature to T > A, this effect 
becomes smeared out by thermal broadening of the Fermi 
distribution, and the injection rate decreases as shown in 
the insets to Figs. [2] and [3] 

Similar arguments can be also applied to the transi¬ 
tions between n = 1 and n = 4 subbands. As a result, 
one gets a relatively small negative peak in the injection 
rate at Hco ~ 2y, see Fig. [2] (b). The weakness of this in¬ 
jection channel is due to a relatively small Fermi velocity 
in the subband 4 at /j - a •€ a, while its reversed sign is 
due to the opposite spin orientation in these subbands, 
which results (similar to Eq. (Illall and Fig. 0]) in a dif¬ 
ferent shape of the Fermi surface, where the transitions 
begin to occur at k y < 0. In the limit a <C y, the positive 
and negative contributions compensate each other, and 
the current injection efficiency tends to zero, as expected 
in the absence of spin-orbit coupling. 

Having discussed the strong peaks in the current injec¬ 
tion rate, let us consider now briefly the broad structure. 
It is formed by momentum dependence of the matrix el¬ 
ements and velocity, and has the efficiency of the order 
of A/a. The current injection stops at Hlo ss 2y + a, 
where the contributions due to transitions between differ¬ 
ent bands compensate each other. We also mention that 
for B || x , the charge current has only the y-component 
for both polarizations of the incident light. 

Now let us address briefly the problem of spin and spin 
current injection. For both incident light polarizations 
one obtains a net spin polarization along the x and y axes. 
The numerical results are presented in Fig. [5] (a) for the 
total spin polarization S x , while Fig. 0(b) shows injected 
spin polarization associated with specific optical transi¬ 
tions. Physical mechanism of the optically injected spin 
polarization is rather clear, since the spin-flip transitions 
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Figure 5: (Color online), (a) Total injected normalized spin 
polarization S x = X!n n ' ■ Here the Rashba SO cou¬ 

pling a = 4 meV (solid red line) and a = 13 meV (dashed 
red line), y = 5 meV and B = 5 T. The orientations of A 
and B are parallel to the x— axis, (b) Transition-related spin 
injection . 

are related to the above-mentioned fact that the eigen¬ 
states of Hamiltonian (ID are not the spin eigenstates, 
and the broken in magnetic field time-reversal symme¬ 
try allows one to inject spin density. Since the charge 
current is along the y-axis, we obtain effectively a spin- 
polarized current transferring in-plane spin components 
in the y— direction. As concerns the spin current defined 
in Eq. m, it is symmetric with respect to the time re¬ 
versal and, therefore, magnetic field produces there only 
changes proportional to B 2 . 

Summary: We have calculated optical injection of 
charge current in graphene as the photogalvanic effect 
due to spin-orbit coupling [25^ . The current is injected 
only in a finite range of infrared light frequencies, de¬ 
termined by the chemical potential y and the spin-orbit 
coupling strength. The striking feature of the injection 
is a narrow peak at the resonant frequency tux ss 2y, 
where the current injection can be very efficient. Com¬ 
paring the w-dependence of the current and spin injec¬ 
tion, we conclude that, depending on the light frequency, 
one can inject either spin-polarized net electric current or 
net spin polarization without the current injection. This 
result can be applied to a controllable current generation 
in spin-orbit coupled graphene. 
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